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A Monte Carlo study of the polar and azimuthal angular distributions of the 
lepton pair arising from the decay a W or Z boson produced at high transverse 
momentum in hadronic collisions is presented. In the absence of cuts on the final 
state leptons, the lepton angular distribution in the gauge boson rest frame is de- 
termined by the gauge boson polarization. Numerical results for the lepton angular 
distributions in the Collins-Soper frame with acceptance cuts and energy resolution 
smearing applied to the leptons are presented. In the presence of cuts, the lepton 
angular distributions are dominated by kinematic effects rather than polarization 
effects, however, some polarization effects are still observable on top of the kinematic 
effects. Polarization effects are highlighted when the experimental distributions are 
divided by the Monte Carlo distributions obtained using isotropic gauge boson decay. 



I. INTRODUCTION 



The measurement of the angular distribution of leptons from the decay of a gauge boson 

V I1I2 \y = W, Z] produced in hadronic colhsions via a Drell-Yan type process hi + h2 — >■ 

V + X provides a detailed test of the production mechanism of the gauge boson. In the 
absence of cuts on the final state leptons, the lepton angular distribution in the gauge boson 
rest frame is determined by the gauge boson polarization. At 0{as) in perturbative QCD 
the angular distribution is described by six helicity cross sections, which are functions of the 
transverse momentum and rapidity of the gauge boson. The lepton angular distribution can 
be used to analyze the gauge boson polarization and thereby test the underlying production 
dynamics in much more detail than is possible by rate measurements alone. 

The hadronic production of a vector boson V is described in lowest order [0(a°)] by the 
Drell-Yan quark-antiquark annihilation process qq — * V. In the case of W boson production 
at low transverse momentum pj,, the W bosons are produced fully polarized along the beam 
direction due to the pure V — A coupling of the charged currents, and the resulting angular 
distribution of the charged lepton in the W boson rest frame is simply 

dN/d{cos9) ~ valence quarks ® (1 + cos 6')^ + sea quarks ® (1 + cos 6'^) . (1) 

Here 6 denotes the angle between the electron (positron) and the proton (antiproton) di- 
rection. The measurement of this angular distribution^ for the W boson by UAl [||] and 
UA2 1^ established the V — A coupling and the spin-1 nature of the W boson. For Z boson 
production at small prp, the angular distribution can be parametrized as 

dN/d{cos ^) ~ 1 + cos ^ + cos^ 9 . (2) 

The coefficient A^ depends on the vector and axial-vector couplings of the fermions to the Z 
boson and is thus sensitive to sin^ 9w- Unfortunately, measurements of the cos 9 distributions 



At the CERN SppS collider energies the contribution from sea quark annihilation is negligible. 
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for gauge bosons produced at low are only sensitive to the trace of the density matrix 
elements of the vector bosons. 

For gauge bosons produced with transverse momentum [balanced by additional gluons 
or quarks], the event plane spanned by the beam and gauge boson momentum directions 
provides a convenient reference plane for studying the angular distributions of the decay 
leptons. The angular distributions are now sensitive to non-diagonal elements of the gauge 
boson density matrix elements. In leading order QCD [0(0;^)] the angular distribution has 
the general form 

a ~ (1 + COS02) + -Ao(l -Scos^^) + Ai sin 26 cos (h 

+ - sin^ 9 cos 2(j) + ^3 sin 9 cos (j) + cos 9 , (3) 

where 9 and (f) denote the polar and azimuthal angle of the decay leptons in the gauge boson 
rest frame. Leading order (LO) analytical results for the coefficients Ai for the decay of a ly 
or Z boson in the Collins-Soper frame are discussed in Refs. The complete next- 

to-leading-order (NLO) corrections to the parity conserving coefficients were calculated in 
Ref. and were found to be fairly small [less than 10%]. This is because the coefficients Ai 
are ratios of helicity cross sections [see Eq. (^)] and the QCD corrections tend to cancel in 
these ratios. LO results for the coefficients Aq, Ai, and A2 for a virtual photon^ 7* decaying 
into two leptons are discussed in Ref. ||^. An investigation of the decay lepton spectrum in 
the laboratory frame can be found in Ref. [Qj. 

In this paper we present a numerical study of the decay lepton angular distribution 
for W and Z boson production at the Fermilab Tevatron collider [\fS =1.8 TeV]. The 
calculation includes acceptance cuts on the decay leptons. These cuts are necessary to 
reduce the background from low p-p jets which can fake electrons. They are also necessary 
due to the finite acceptance of the detector. Uncertainties due to the finite energy resolution 
of the detector are simulated in the calculation by Gaussian smearing of the lepton four- 



^The parity- violating coefficients ^3 and vanish in this case. 
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momentum vectors. The energy resolution smearing has a non-neghgible effect only on the 
cos^ distribution of the lepton from W decay. We show the transverse momentum [Pt(V")] 
dependence of the coefficients Ai and study the effect of the acceptance cuts and energy 
resolution smearing on the (p aiid cos 6' distributions in the Collins-Soper (CS) frame. The 
CS frame has the unique advantage that the polar and azimuthal angles 6 and can be 
reconstructed modulo a sign ambiguity in cos^^ without information on the longitudinal 
momentum of the neutrino. 

Since the effects of the NLO corrections are small, it is sufficient to use LO matrix 
elements in our Monte Carlo study. At NLO there are three more angular coefficients in 
Eq. which receive contributions from the absorptive part of the one-loop amplitudes. A 
numerical analysis shows that their contribution is small [at the level of 1% at high pj] 0, 
thus they will be neglected in this analysis. 

The measurements of the lepton distributions arising from the decay of a virtual photon 
7* in nucleon-nucleon scattering presented in Ref. P appear to disagree with the QCD im- 
proved parton model. Several different nonperturbative effects have recently been discussed 



in Ref. [|TU| to explain the discrepancy. It would be interesting to test whether the angular 
distributions from W and Z boson decay at the Tevatron are in agreement with theoretical 
predictions. 

The remainder of this paper is organized as follows: the formalism for describing the 
angular distributions is discussed in Sec. II, numerical results are presented in Sec. Ill, and 
summary remarks are given in Sec. IV. 

II. ANGULAR DISTRIBUTIONS 

We consider the angular distribution of the leptons coming from the leptonic decay 
of gauge bosons produced with non-zero transverse momentum in high energy proton- 
antiproton collisions. For definiteness we take 

p(Pi) + p{P2) ^W^{Q)+X^ l^ih) + uiik) + X , (4) 



and 

p(Pi) + p{P2) ^ Z{Q) + X ^ + /- {h) + X , (5) 

where the quantities in parentheses denote the four-momenta of the particles. In leading 
order QCD [C(a;s)] the parton subprocesses 

q + q-^V + g, 

q + g-^V + q, (6) 
g + q^V + q, 

contribute to high pj, gauge boson production. 

In the parton model the hadronic cross section is obtained by folding the hard parton 
level cross section with the respective parton densities: 

where the sum is over a,b = q, q, g. (x, M^) is the probability density for finding parton 
a with momentum fraction x in hadron h when it is probed at scale fi'jp. The parton level 
cross section for the processes in Eq. (|^) are denoted by d&ab- Denoting hadron level and 
parton level quantities by upper and lower case characters, respectively, the hadron and 
parton level Mandelstam variables are defined by 

S={P, + P2)\ T={P,-Q)\ U = {P2-Q)\ (8) 

and 

s = +P2y = X1X2S , 

t = {Pi-Q? = x^{T-Q^)+Q\ (9) 
u = {P2-Q? = X2{U-Q')+Q\ 

where pi = xiPi and p2 = xiP2- The rapidity y of the gauge boson in the laboratory frame 
can be written 
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and the transverse momentum prp of the gauge boson is related to the Mandelstam variables 
via 

= Wl^i^iWl^ - = . (11) 

5* s 

The angles 6 and (p in dQ* = d cos 6 d(f) are the polar and azimuthal decay angles of the 
leptons in the gauge boson rest frame measured with respect to a coordinate system to be 
described later. 

Technically, the lepton-hadron correlations are described by the contraction of the lepton 
tensor L^^, with the hadron tensor H^'^, where L^,^ acts as an analyzer of the gauge boson 
polarization. The angular dependence in Eq. (|^ can be extracted by introducing helicity 
cross sections corresponding to the non-zero combinations of the polarization density matrix 
elements 

H^„,, = e;{m) e,{m') , (12) 

where m, m' = +1, 0, —1 and 

e,{±l) =^(0;±l,-z,0), 

(13) 

e^(0) =(0; 0,0,1), 

are the polarization vectors for the gauge boson defined with respect to the chosen gauge 
boson rest frame. The angular dependance of the differential cross section can be written 
[see Ref. for details] 

IGtt da da^^^ , ^ ^. da^ , , 

1 + cos^g + , o , l-3cos^^ 



3 dj)\ dy d cos 9 d(j) dpj. dy dpj. dy 

dcF _ . 2 n , da^ 
, ^ , 2 sm d cos 20 + —5- 
dp^ dy dp^ dy 



+ ^ ^ ^ 2 sin^ e cos sin26'cos<^ (14) 



d(T^ d(T^ 
+ „ 4\/2 sin 6 cos + „ 2 cos 9 , 
dp J' dy dp J. dy 



where 



a 



(15) 



5 



The unpolarized differential production cross section is denoted by cr^"^'^ whereas a^'^'^'^'^ 
characterize the polarization of the gauge boson, e.g., the cross section for longitudinally 
polarized gauge bosons is denoted by a^, the transverse- longitudinal interference cross sec- 
tion by 0"^, the transverse interference cross section by cr^, etc (all with respect to the 2;-axis 
of the chosen lepton pair rest frame; for details see appendix C of Ref. 0). At NLO there 
are three more "T-odd" angular coefficients in Eq. (|T¥) for W and Z production. However, 



their numerical contribution is small |^,n| and we will neglect them in the present analysis. 
The helicity cross sections a" contain the following coupling coefficients: 

^^+"'™~(^? + anK + «?), (16) 
cr-^'^ ~ vi ai Vq ag , (17) 

where Vg{vi) and ag{ai) denote the vector and axial vector coupling of the gauge boson to 
the quark (lepton). 

The hadronic helicity cross sections in Eq. (JT^) are obtained by convoluting the 

partonic helicity cross sections with the parton densities: 

^ = / rfx,rfx./'^Hx„;.|)/'^^(x„/.|)^. (18) 

The helicity cross sections a" are dependent on the choice of the z-axis in the rest frame of 



the gauge boson and are explicitly given in Ref. for the Collins-Soper frame. In this 
frame the z-axis bisects the angle between Pi and —P2'- 



Pi = El (sin 7^5, 0, cos 7^5) , 
P2 = E2 (sin 7^5, 0, - cos 7^.5) 



(19) 



with 



cos -ics 



1/2 



1/2 



XT-Q^)iU-Q^) 



El 



Eo 



(20) 

(21) 
(22) 



In the laboratory frame, the 2;-direction is defined by the proton momentum and the x- 
direction is defined by the transverse momentum of the gauge boson. The Lorentz transfor- 
mation matrix for the boost from the laboratory frame to the CS frame is given by 



Px 

Py 



Pz 

V J 



( Qo 



/Q2 



Pt Qo 







Qz 











1 



^Xt 





Q Qo 

Xt 





(e\ 




Px 




Py 




Pz 



(23) 



CS 



V / 



lab 



where 



Q^+pI 



(24) 



Concerning the leptonic decay modes, W bosons have the disadvantage that the decay 
kinematics cannot be completely reconstructed due to the unobservability of the neutrino. 
The transverse components of the neutrino's momentum can be approximated from trans- 
verse momentum balancing in calorimetric experiments, however, the longitudinal compo- 
nent can only be obtained up to a twofold ambiguity. By choosing the Collins-Soper frame 
0, the polar angle {0) and the azimuthal angle (0) of the charged lepton can be reconstructed 
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modulo a sign ambiguity in cos 6, thereby enabling the measurement of four [cr'^+^'^'^'^] out 
of the six helicity cross sections in Eq. (p!^ |0. For Z boson production and decay, (j) ^^'^ 
cos 9 can be determined without ambiguity. 

Let us briefly review how 9 and cj) can be obtained from measured quantities in W boson 
production and decay. Denote the charged lepton's four-momentum in the laboratory frame 
by Vi = {Ei,pix,piy,piz) where pi^ and pix are the lepton's momentum components parallel 
and perpendicular to the proton direction in the event plane, i.e., the plane spanned by the 
proton and the W boson with pj. > 0. The lepton's four-momentum in the CS frame, pp'^, 
can be obtained by Lorentz-boosting pi from the laboratory frame to the CS frame. One 
obtains 



P?y = Ply , (25) 



Plz 



where the ly-mass constraint has been imposed on the lepton-neutrino system. The ± signs 
correspond to the two solutions for the longitudinal momentum of the neutrino. The trans- 
verse component pf^^ is uniquely determined by the measurable laboratory frame quantities^ 
pix and pj^ = pix + Pux- The CS frame is thus unique in the sense that in this frame the 
lepton's transverse momentum is independent of the unmeasured longitudinal momentum 
of the neutrino. On the other hand note that the longitudinal component pp'^ is determined 
only up to a sign. The angles 9 and are obtained from the charged lepton's momentum 
components in Eq. (|25|) via = arctan (pfy^ /Ptf^ and 9 = arctan ^Jip'tJ'Y + {p'ty^Y / Ph^ ■ 
The twofold ambiguity in the reconstruction of the lepton's longitudinal momenta in the CS 
frame translates into an ambiguity ^ vr — ^ in the polar angle, while is completely de- 
termined. This implies that only the helicity cross sections and in Eq. (p!4D 



Note that there is a misprint in the definition of pj. in Ref. 
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can be determined in a boson production experiment. For Z boson production, all six of 
the angular coefficients in Eq. (|1^ can be determined. 
Introducing the standard angular coefficients 

° ~ d^ ' ' ~ da^+L ' ' ~ d^ ' 

A^da^ ^ 2daP 

the angular distribution in Eq. (14) can be conveniently written 



da 3 d(y^+^ 



dp^dy d cos 9 d(j) IGvr dp^dy 



;i + cos2 0)+ ^Ao (l-3cos^ 



+ Ai sin 2^ cos + -A2 sin^6'cos20 

+ ^3 sin 6^ cos + COS 6^ ] . (27) 

Integrating the angular distribution in Eq. ( pTj) over the azimuthal angle yields 

da 

dp\ dy d cos 6 

where 



C(l + aicose + a2Cos^^) , (28) 



C = - 



3 da 



U+L 



8 dpj' dy 

Integrating Eq. ([27| ) over the polar angle 9 yields 

da 1 da^+^ 



2^4 2-3Ao 



dp^dydcj) 271 dp^dy 

where 



1 + /^i cos + /^a cos 20) , (30) 



ft = |43, A = ^. (31) 

Before discussing numerical results, let us briefly discuss a possible strategy for extract- 
ing the angular coefficients. By taking moments with respect to an appropriate product 
of trigonometric functions it is possible to disentangle the coefficients Ai. A convenient 
definition of the moments is 
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^ / dajpT, y, 9,(j)) md cos 9 d(j) 
J da {p y, 9, (f)) d cos 9 d(j) 



which leads to the following results: 



;i) = 1 , (33) 



;i(l-3co.= »)) = l (Ao-|). (34) 

(sin2^cos0) = - , (35) 

(sin2^cos20) = ^ ^2, (36) 

(sin^cos0) = i A3 , (37) 

(cos^) = ^A4. (38) 



III. NUMERICAL RESULTS 

In this section numerical results are presented for high pj. production and leptonic decay 
of W and Z bosons at the Tevatron collider center of mass energy [v^ = 1.8 TeV]. The nu- 
merical results have been obtained using the MRS set D'_ parton distribution functions 



with Aj^ = 215 MeV and the one- loop formula for ag. The renormalization scale /i^ and 
the factorization scale fi'jp in Eq. (^ have been taken to be /^/j = /i^ = My + p'^{V), where 
My and Pj-iV) are the mass and transverse momentum, respectively, of the gauge boson. 

We begin with numerical results for the coefficients Ai in Eq. (pTl). These coefficients have 
been extracted from the Monte Carlo program by using the moments defined in Eq. (p2D. 
Figure 1(a) shows the coefficients Aq, A2, and ^3 in the CS frame as functions of prpiW) for 
W production. The angles 9 and (p are the decay angles of the charged lepton in the CS 
frame. The results are identical for and W~ . The lepton's four-momentum in the CS 
frame was obtained using Eq. ([25|). Note that because of the ambiguity 9 ^ n ~ 9 in the 
polar angle reconstruction in the CS frame, the coefficients Ai and are not observable 
in a boson production experiment. However, all of the A^ coefficients are observable for 
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Z boson production, and the numerical results are shown in Fig. 1(b), where Q and are 
the angles of the negatively charged lepton. No acceptance cuts have been applied to the 
leptons. 

The coefficients and A2 are increasing functions of Vt^^ ^^"^ ^^e deviations from the 
lowest order expectation of Eq. (0) [Aq = = 0] are quite large, even at modest values of 
■pj.{y\ i.e., for Py(V^) ~ 20 — 50 GeV. It has been noted in Ref. that these coefficients 
are exactly equal in LO. This is no longer true in NLO, but the corrections to Aq and A2 are 
fairly small It has been shown in Ref. that the 0{as) relation Aq = A2 is peculiar to 



a vector gluon theory. For a scalar gluon this relation is badly broken, e.g., Aq — A2 ^ —2 
in the Gottfried- Jackson frame at the CERN SppS energies. In was in fact the study of this 
relation in W decays by UAl which fixed the spin of the gluon [|l|. NLO numerical results 



for Aq and A2 for W decay in the Gottfried- Jackson frame are given in Ref. [|14| and the 
deviations from Aq = A2 are found to be small. NLO numerical results for Aq and A2 in the 
CS frame are presented in Refs. [|6|,p!5|. 

The deviation of Ai and A^ from zero is much smaller, even at large pj.{V). This is 
a special effect of the CS frame. In this frame the qq contribution in Eq. (j^) to Ai and 
A3 is antisymmetric under the interchange of Xi and X2', this can be explicitly seen in the 
corresponding matrix elements in Eq. (24) of Ref. The qg process gives a sizable con- 
tribution to ^3 only for W boson production, and a measurement of ^3 could in principle 
provide constraints on the gluon distribution function^. However, as we will see later, exper- 
imental cuts introduce additional complicated angular effects and the resulting data sample 
can no longer be described by the simple angular distribution in Eq. (p?)). Futhermore, the 
extraction of the different polarization cross sections via the moments in Eq. (|32|) becomes 
problematic. Note that A4^ in Fig. 1(b) is proportional to vf af and is therefore sen- 

sitive to the Weinberg angle 6^. This determination would be totally independent of the 



^Note that A-^ has the wrong sign in the figures in Ref. 
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W and Z boson masses. However, a precise determination of sin^ 9w by this method is also 
hmited by the problem of extracting with sufficient accuracy when experimental cuts are 
applied to the decay leptons [see the discussion below]. 

To give a feeling for the and cos^ distributions of the decay leptons in the CS frame 
we show numerical results for the coefficients ai,a2 and /3i,/?2 [see Eqs. ( PB] ) and (^)] as a 
function of prp{V) for the W and Z boson in Fig. 2. As in Fig. 1, no cuts have been applied 
to the leptons. The coefficients are again very sensitive to the transverse momentum of 
the gauge boson. In the remainder of this paper we investigate the feasibility of using this 
dependence to test the polarization properties of the W and Z bosons at the Tevatron. 

We want to point out that we are not in agreement with the theoretical predictions for the 
angular coefficients of the Z boson presented in Ref. [Q where the "soft gluon resummation" 
formalism was used for the integrated cross section^ but not for the individual helicity 
cross sections which are responsible for the angular distributions. Even for gauge boson 
transverse momenta larger than pj, > 20 GeV, where the angular coefficients should be 
reliably predicted by fixed order perturbative QCD, our predictions differ dramatically from 
the results in Ref. [^, for example, by a factor five for the coefficient Aq in Eq. ( [27| ) at 
p-riZ) = 60 GeV. The complete NLO calculation in Ref. [^] for the helicity cross sections 
shows that the deviations from the LO result are much smaller than the predictions in 
Ref. §. 

Figures 3-13 show the normalized and cos 6' distributions for the leptons from the 
decay of W and Z bosons for four bins in the transverse momentum of the gauge boson. 
To demonstrate the effects of acceptance cuts, results are shown first without cuts and then 
with typical acceptance cuts imposed on the leptons. The acceptance cuts are necessary to 
reduce the background from low pj, jets which can fake electrons. They are also needed due 



^The integrated cross section is essentially the total cross section after integration over 9 and 0, 
i.e., a^+^ in Eq. (ITl). 
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to the finite acceptance of tlie detector. Measurements of the charged lepton momentum and 
the missing transverse momentum have inherent uncertainties due to the finite energy 
resolution of the detector. These uncertainties have been simulated in our calculation by 
Gaussian smearing of the charged lepton and neutrino four-momentum vectors with standard 
deviation a. The numerical results presented here were made using a values based on the 
CDF specifications [|l^. The energy resolution smearing has a non- negligible effect only on 



the cos distribution of the charged lepton from W decay. 

Figure 3 shows the and cos 9 distributions for the W boson without cuts or smearing. 
The curves in Fig. 3 can be obtained from a2, and (32 in Fig. 2. For example, in the 
lowest Pxi^) bin [solid curve], the distribution is almost flat [/3i,/32 ~ 0] and the cos^ 
distribution is approximately 1 + cos 6*^ [^2 ~ !]• In the highest priW) bin [dot-dashed 
curve] , one observes the dependence [/3i , 7^ 0] resulting from the non-diagonal elements 
of the spin-density matrix of the W boson, whereas the corresponding cos 6' distribution is 
almost fiat [0:2 ~ 0]. Note that if the W boson was to decay isotropically, the and cos^ 
distributions would both be fiat. 

Figure 4 shows the and cos^ distributions for the W boson for the same bins in prp{W) 
as in Fig. 3, but with energy resolution smearing and the cuts 

Pt{1) > 25 GeV, > 25 GeV, y{l) < 1. (39) 



Only 33% of the events pass the cuts in Eq. (pOj). The cuts have a dramatic effect on 
the shapes of the distributions. The shapes of the distributions are now governed by the 
kinematics of the surviving events. The cuts, which are applied in the laboratory frame, 
introduce a strong dependence. The "kinematical" distribution in Fig. 4(a) is very 
different from the "dynamical" distribution in Fig. 3(a). The only remaining vestiges of 
the polarization effects in the distribution are the dips in the high Px{W) curve [dot-dashed 
curve] at = 90°, 180°, and 270°. 

The cos 6' distributions with cuts in Fig. 4(b) are very different from the corresponding 
results without cuts or smearing in Fig. 3(b). In order to separate the effects of smearing 
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from the effects of cuts, the cos^ distribution is shown in Fig. 5 with the cuts of Eq. (p9D, 
but without smearing. The main effect of the cuts is to remove events with cos6 > 0.5 
for small values of p^iW), causing the curves to drop to zero as cos 6* — > 1 [see the solid 
and dashed curves in Fig. 5]. Comparing Figs. 4(b) and 5, we see that the smearing of 
the lepton momentum and the missing transverse momentum has a large effect on the cos 6 
distribution in the small cos^ region. This behavior of the smeared curves can be traced 
to the expression for pf/ in Eq. (|25|) . Small values of cos^ correspond to small values of 
pf/. When events with small values of cos^ are smeared, the argument of the square root 
function in the expression for p^"^ often becomes negative, in which case the argument is 
set equal to zero in the Monte Carlo program. Thus many events with small values of cos 6 
are smeared such that they end up in the cos 6* = bin, which is not shown here since it 
is off the vertical scale. When cuts and smearing are combined, as in Fig. 4(b), very little 
polarization dependence is left in the cos 6 distribution for the W boson. 

We have also analyzed the effect of the cuts separately by using the correct matrix 
element for W boson production, but with isotropic decay of the W boson, i.e., neglecting 
spin correlations between W production and decay. The angular distributions in this case are 
very similar to the ones shown in Fig. 4 for the full matrix element; the remnant polarization 
effects discussed in Fig. 4 are of course absent. In Fig. 6 we show ratios of the and 
cos 6' distributions for the same bins in prp{W) as in Figs. 3-5; the distribution with full 
polarization has been divided by the distribution obtained with isotropic decay of the W 
boson. Cuts and smearing are included in both cases. The large effects from the cuts 
and smearing are expected to almost cancel in this ratio. In fact, we nearly recover the 
and cos 6' dependence of Fig. 3 which contains no cuts or smearing. Therefore, to regain 
sensitivity to the polarization effects in the presence of large kinematic cuts, we propose 
to divide the experimental distributions by the Monte Carlo distributions obtained using 
isotropic gauge boson decay. 

Figure 7 again shows the and cos 6 distributions for the W boson with energy resolution 
smearing, but now with the looser cuts 
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Pt{1) > 15 GeV, ^j, > 15 GeV, y{l) < 2.5 . (40) 

The results now display more of the polarization effects seen in Fig. 3. The dips in the (j) 
distribution in the high priW) curve [dot-dashed curve] at = 90° and 270° in Fig. 7(a) 
are much more distinctive than in Fig. 4(a). Dynamical effects in the cos^ distribution [see 
Fig. 3(b)] can still be observed in Fig. 7(b) for cosO values in the range 0.3 < cos^^ < 0.8. 
Figure 8 shows the ratio of the "polarized" (p and cos 6* distributions to the corresponding 
distributions with isotropic leptonic decay for the cuts in Eq. (|^). Most of the polarization 
dependence seen in Fig. 3 is retained in this ratio. 

The and cos^ distributions of the leptons from Z boson decay are shown in Figs. 9-13 
for four bins in prp[Z). For Z bosons, the lepton momenta in the CS frame can be obtained 
from the measured lepton momenta in the laboratory frame by applying the boost matrix 
given in Eq. (^3]). The (p and cos^ distributions of the negatively charged lepton are shown 
without cuts in Fig. 9. Again, there are large differences between the different Px{Z) bins. 
If the electric charge of the lepton can be determined, then the coefficient A4 [or ai] can in 
principle be measured and the cos^ distribution in Fig. 9(b) is asymmetric about cos6 = 0. 
However, if the lepton's electric charge can not be determined, the term linear in cos^ will 
be averaged out, and the resulting cos 6* distribution will be symmetric. 

Figure 10 shows the (p and cos^^ distributions for the Z boson using the same bins in 
Pt{Z) as in Fig. 9, but now with energy resolution smearing and the cuts 

p^il) > 25 GeV, y{l) < 1. (41) 

Note that these cuts are more stringent than the cuts for the W boson in Eq. ([39| ) since 
now both leptons must be in the central rapidity region. The effect of the energy resolution 
smearing is negligible for the Z boson. The cuts introduce a (p dependence similar to that 
discussed for the W boson case in Fig. 4(a). Some polarization effects are still observable 
in the (p distribution in the high Pt{Z) bins, for example, the dip in the dot-dashed curve 
and the flat behavior in the dotted curve around cp = 90° and 270° in Fig. 10(a) are due to 
polarization effects. 
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The cuts in Eq. (^) have an even more dramatic effect on the cos 9 distribution. Since 
nearly all events with | cos9\ > 0.6 are rejected, almost all sensitivity to the asymmetry in 
Fig. 9(b) is lost. Futhermore, the curves in Fig. 10(b) are fairly similar in shape. Therefore 
it may be very difficult to observe any polarization effect from the Z boson in these distri- 
butions when the cuts in Eq. (^Tj) are imposed. However, by forming the ratio of the and 
cos^ distributions with full polarization to the corresponding distributions obtained with 
isotropic Z decay, most of the polarization dependence is recovered; see Fig. 11. 

Since the experimental signature of the Z boson decay is much cleaner than that for the 
W boson. Fig. 12 shows the and cos 6 distributions for the Z boson for the looser cuts 

Pt{1) > 10 GeV, y{l) < 2.5 . (42) 

The large dips in the distribution in the high p^iZ) bin [dot-dashed curve] at = 90° 
and 270° are due to polarization effects [compare Fig. 12(a) to Fig. 9(a)]. The curves in 
the central region of cos^ in Fig. 12(b) now clearly exhibit the polarization effects seen in 
the corresponding curves in Fig. 9(b). Finally, Fig. 13 shows the ratio of the "polarized" 
and cos^ distributions to the corresponding distributions obtained with isotropic leptonic 
decay for the Z boson for the cuts in Eq. (|4^) . The ratios once again contain most of the 
polarization dependence seen in Fig. 9. 

IV. SUMMARY 

The polar and azimuthal angular distributions of the lepton pair arising from the decay 
of aW or Z boson produced at high transverse momentum in hadronic collisions have been 
discussed. In the absence of cuts on the final state leptons, the general structure of the 
lepton angular distribution in the gauge boson rest frame is determined by the gauge boson 
polarization. At 0{as) in perturbative QCD, the structure is described by six helicity cross 
sections, which are functions of the transverse momentum and rapidity of the gauge boson. 

We have studied the angular distributions of the leptonic decay products of high pj, 
gauge bosons when acceptance cuts and energy resolution smearing are applied to the lep- 
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tons. When acceptance cuts are imposed on the leptons, the shapes of the lepton angular 
distributions are dominated by kinematic effects and the residual dynamical effects from 
the gauge boson polarization are small. The kinematic effects become more dominate as 
the cuts become more stringent. Energy resolution smearing has a significant effect only 
on the cos^^ distribution from W decay. The large smearing effect in this distribution is 
a consequence of the undetermined longitudinal momentum of the neutrino. The angular 
distributions have been calculated in the Collins-Soper frame. For W decay, this frame has 
the unique advantage that the polar and azimuthal angles 6 and can be reconstructed 
modulo a sign ambiguity in cos 6 without information on the longitudinal momentum of the 
neutrino. 

When cuts are imposed on the leptons, the only polarization effects visible in the 
and cos^ distributions occur at large values of the gauge boson transverse momentum, 
where unfortunately, the cross section is smallest. Polarization effects can be maximized by 
minimizing the cuts, however, this strategy is severly limited since cuts are needed to reject 
background events. Alternatively, it may be possible to highlight gauge boson polarization 
effects by "dividing out" the kinematic effects, i.e., if the histogrammed data is divided by 
the theoretical result for isotropic gauge boson decay, the resulting ratio is more sensitive 
to polarization effects. 

Since polarization effects are nearly obscured by the kinematical effect of cuts, experi- 
mental analyses of the very small so-called "T-odd" effects discussed in Ref. or possible 
CP violation effects in the Standard Model through the Kobayashi-Maskawa mechanism in 
hadronic W and Z production [|1^ appear to be impractical in the presence of realistic cuts 



at the Tevatron. 
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FIGURE CAPTIONS 



Fig. 1 a) Angular coefficients Aq, A2, and A3 for W boson production and decay in the CS 
frame as a function of the transverse momentum prp{W) at =1.8 TeV. 
b) Angular coefficients Ao,Ai,A2,A3, and A4 for Z boson production and decay in 
the CS frame as a function of the transverse momentum Pt{Z) at y/S — 1.8 TeV. 
No cuts or smearing have been applied. 

Fig. 2 a) Angular coefficients q;2,/3i, and P2 for W boson production and decay in the CS 
frame as a function of the transverse momentum P'j'{W) at =1.8 TeV. 
b) Angular coefficients ai,a2, Pi, and (^2 for Z boson production and decay in the CS 
frame as a function of the transverse momentum Pt{Z) at = 1.8 TeV. 
No cuts or smearing have been applied. 

Fig. 3 a) Normalized <p and b) normalized cos^ distributions of the charged lepton from W 
boson decay in the CS frame. Results are shown for four bins in Pt{W): 
10 GeV < Pt{W) < 20 GeV (solid), 
20 GeV < p^iW) < 30 GeV (dashed), 
30 GeV < p^iW) < 70 GeV (dots), 
70 GeV < Pt(W) (dot-dashed). 
No cuts or smearing have been applied. 

Fig. 4 Same as Fig. 3 but with smearing and the cuts p^il) > 25 GeV, > 25 GeV, and 
y{l) < 1. 

Fig. 5 Normalized cos^^ distribution of the charged lepton from W boson decay in the CS 
frame. The cuts Pt{1) > 25 GeV, > 25 GeV, and y{l) < 1 have been imposed, but 
no smearing is included. 

Fig. 6 Ratios of distributions obtained with full polarization effects to those obtained with 
isotropic decay of the W boson. Parts a) and b) are the ratios for the ^ and cos^ 
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distributions, respectively. Energy resolution smearing and the cuts Pr(0 > 25 GeV, 
> 25 GeV, and y{l) < 1 are included. 

Fig. 7 Same as Fig. 3 but with smearing and the cuts Pt{1) > 15 GeV, > 15 GeV, and 
y{l) < 2.5. 

Fig. 8 Same as Fig. 6 but with the cuts > 15 GeV, > 15 GeV, and y{l) < 2.5. 

Fig. 9 a) Normahzed and b) normalized cos 9 distributions of the negatively charged lepton 
from Z boson decay in the CS frame. Results are shown for four bins in Pj-{Z): 
10 GeV < Pt{Z) < 20 GeV (sohd), 
20 GeV < Pt,{Z) < 30 GeV (dashed), 
30 GeV < Pt{Z) < 70 GeV (dots), 
70 GeV < Pt{Z) (dot-dashed). 
No cuts or smearing have been applied. 

Fig. 10 Same as Fig. 9 but with smearing and the cuts Pt{1) > 25 GeV and y{l) < 1. 

Fig. 11 Ratios of distributions obtained with full polarization effects to those obtained with 
isotropic decay of the Z boson. Parts a) and b) are the ratios for the (j) and cos 9 
distributions, respectively. Energy resolution smearing and the cuts Pr(/) > 25 GeV 
and y{l) < 1 are included. 

Fig. 12 Same as Fig. 9 but with smearing and the cuts Pt{1) > 10 GeV and y{l) < 2.5. 

Fig. 13 Same as Fig. 11 but with the cuts prp{l) > 10 GeV and y{l) < 2.5. 
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